Butterfly-like spectra and collective modes of antidot superlattices in magnetic fields 
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We calculate the energy band structure for electrons in an external periodic potential combined 
with a perpendicular magnetic field. Electron-electron interactions are included within a Hartree 
approximation. The calculated energy spectra display a considerable degree of self-similarity, just 
as the "Hofstadter butterfly." However, screening affects the butterfly, most importantly the band- 
widths oscillate with magnetic field in a characteristic way. We also investigate the dynamic response 
of the electron system in the far-infrared (FIR) regime. Some of the peaks in the FIR absorption 
spectra can be interpreted mainly in semiclassical terms, while others originate from inter(sub)band 
transitions. 
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Recent years have witnessed a considerable amount 
of research effort directed towards understanding of 
the physics of two-dimensional electron systems (2DES) 
whose dimensionality is further restricted by man-made 
periodic potentials and perpendicular magnetic fields. 
These include quantum dot arrays and antidot superlat- 
tices. Concentrating on the latter ones, one distinguishes 
two principal directions of experimental work: transport 
studies and far infrared (FIED spectroscopy. Some of 
the transport measurementsBH have been performed in 
search of evidence for a self-similar energy spectrum, the 
so called Hofstadter butterflyE A main theme in the FIR 
absorption experiments has been to detect and classify 
the richj-variety of collective modes that occur in these 
systems dQ 

Along with the experigipital work, theorists have ad- 
dressed the same issues. 13 The main difficulty lies in 
the fact that while the superlattice is periodic, the Hamil- 
tonian (including a vector potential) is not. Most recent 
calculations of superlattice electronic structure have used 
the Ferrari basis to deal with this matter .OrWfi-will in- 
stead apply ray-group-theoretical techniquesEjTS to ef- 
fectively reduce the calculational complexity. 

With this approach, we are able to go beyond earlier 
calculations. We find the band structure for interact- 
ing electrons in general "rational" magnetic fields [i.e., 
the flux through a unit cell is (L/iV)$o, where $0 is a 
flux quantum and L,N G Z\. Consequently, we are able 
to trace even fine-scale features of the butterfly and at 
the same time study the effects screening has on it. We 
also explore the FIR response. The resulting spectra are 
rather rich. Along with absorption peaks caused by col- 
lective modes, and known from experimentsjju we find 
additional ones of mostly quantum-mechanical origin. 

The antidot superlattice considered here is of simple 
square symmetry R = niai -I- n2SL2, with lattice param- 
eter a. The effective one-particle Hamiltonian is 

if = fp+-Af + V^;(G)e'G■^ (I) 
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gauge) describes the perpendicular magnetic field B and 
G — 5ibi + §2^2, denotes the reciprocal lattice vec- 
tors. We use GaAs parameters and work with short- 
period superlattices with a = 1000 A and electron den- 
sity Us = 1.2 • 10^^ cm~^. A typical magnetic field 
B — 1.65 T gives four flux quanta per unit cell, filling 
factor V — 1.5, and the cyclotron energy TilOc = 2.86 meV, 
where Wc — eB/mc. The last term in Eq. (|l|) is a sum 
of the external superlattice potential, described by a few 
principal Fourier components, Ej and the Hartree poten- 
tial. As for electron spin, we keep the twofold degeneracy 
in mind when counting states, but neglect other effects 
such as Zeeman splitting and exchange interaction. 

The electronic states we set out to solve for will satisfy 
the modified Bloch conditions 



TA/(ai(2))Vq = e-''^i<=)"iVq 



(2) 



when the magnetic flux through a unit cell equals an in- 
teger number of flux quanta $0. Here 



rM(R) = exp 



.1h-(p-£a 



(3) 



are magnetic translation operators forming a ray group, 
and the eigenstates can still be classified by different val- 
ues of the crystal momentum q = gibi -|- q2^2i in the 
first Brillouin zone.E3 Actually, one finds sets of L lin- 
early independent functions each transforming according 
to the same irreducible representation. This manifests 
itself as the splitting of the Landau band into L sub- 
bands. For rational fields with flux {L/N)^f) per unit 
cell, the irreducible representations of the group (||) are 
A^-dimensional and the states iVpfcjldj-degenerate. This 
calls for a generalized treatmentOliZiO which we have 
implemented but do not further describe here. 

The next important step towards a solution is a canoni- 
cal coordinate transformation.llZI We switch to dimension- 
less unitsEj to be used hereafter, and introduce 



where the vector potential A = B x r/2 (symmetric 



S_{-q) = ±Py + x/2, p^(r,) = PxTyf^- 



(4) 
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This preserves the canonical commutators, maps the ki- 
netic energy of the Hamiltonian in the symmetric gauge 
onto a harmonic osciUator in ^, and makes the magneto- 
translations act only on 77, 

Ho^^ + P|) , Tm(R) = exp i-iR^p^ + iRyTj) . (5) 

The periodic potential mixes the ^ and 77 degrees of free- 
dom (in these coordinates it behaves like a magnetic 
translation operator) 

H,=J2^iG)X{G\0Y{G\v), 

G 

X{G\0^eMiG.^-iGyp^), (6) 
Y{G\ri) = exp(iG2;77 + iGyPn). 

Using projection-operator techniques we find the 
symmetry-adapted 7/-dependent functions 

(^(q,;h)= £ e^^J^T^ + ^-^-amV (7) 

m— — 00 ^ ^ 

labeled by the magnetic crystal momentum q in the first 
magnetic Brillouin zone (i.e., < 91,92 < 1) and the 
subband index I = Q, . . . , L — 1. Now the Ansatz 

L — 1 00 

for the eigenstates allows for subband mixing, and the 
^-dependence is accounted for by an expansion in har- 
monic oscillator eigenfunctions Xn- Inserted into the 
Schrodinger equation, the Ansatz yields the eigenvalue 
problem determining the electron states 

+ Y,<G)^n'iG)B„„'iG)}an'i' = Eani. (9) 

G 

The subband and Landau-level mixing coefficients are 

Au>{G) = Q'^-^'^[9i92/2+gil + (qig2~q2gi)+qi(l-l')]/L gimodL) ^ 

/oo 
Xn{OX{G\Oxn'{Odt (10) 
-00 

Equation (^) must be iterated together with the Pois- 
son equation updating the Hartree potential until self- 
consistency is reached. 

Figure |l| shows the splitting of the first four Lan- 
dau levels as a function of the dimensionless inverse flux 
<i>o/(i?a^) = N/L. With a reasonable computational ef- 
fort we could treat rational fields with L < 14 and all 
possible A^'s. This is enough to clearly resolve the intri- 
cate subband clustering.0 



It is easy to see that the bandwidths in Fig. |l| decrease 
with increasing magnetic fiux; however, the decrease is 
not monotonous. Instead they have maxima for fiux val- 
ues 6, 3, and 2, (see the inset) when there are 1, 2, and 3 
completely filled Landau levels, respectively. Then the 
2DES cannot screen the external potential very effec- 
tively, and the Fourier coefficients v{G) are larger than 
for other fiux values. Thus, since the bandwidth is set by 
a competition between the band-narrowing effects of the 
magnetic field and the band broadening tendencies of the 
potential, this leads to a cusped behavior of the band top 
and bottom at integer filling factors. For the filling fac- 
tors v < 1 we also observe the same qualitative behavior 
while quantitative predictions of the Hartree theory in 
this region may be inaccurate. We note that there exist 
other (unrelated to electron-electron interaction) mecha.- 
nisms which also lead to nonmonotonous bandwidths.E^ 

Electron-electron interaction also contributes to dimin- 
ishing the symmetry of the b uttc rfiy as strong coupling 
between different bands does.liiiij The second and third 
bands in Fig. [l] which are traversed by the chemical po- 
tential /i show reduced regularity if compared to well- 
separated noninteracting bands in Fig. 3 (a) of Ref. 
whereas our fourth band, well above /i, would exhibit a 
considerable resemblance to Hofstadter's one-band pic- 
ture when properly rescaled. 

Turning to the dynamic response of the 2DES, we cal- 
culate the density-density response function Rgc (k, uj) 
within the random-phase approximation (RPA) by solv- 
ing the set of equations 

RgG' = PgG' +J2PGG"V,,{k + G")Rg"G'- (11) 
G" 

Here Pgg' is the independent particle response function 
which we can evaluate knowing the electron eigenstates 
already calculated. The FIR absorption of the long wave- 
length (k in the first Brillouin zone) light is proportional 
to -ujlmRoo{u!)E3 

The so calculated spectra typically exhibit several 
conspicuous peaks. In Fig. |^ (a), we display spectra 
calculated for electron densities = 1.2 • 10^^ cm^^ 
and 1.4 • 10^^ cm~^, respectively, and wave vector k — 
{7T/10)a^^x. Following the suggestion of Ref. |l^ to clas- 
sify the different peaks by studying the corresponding 
charge fiuctuations; we also trace their development in 
time to pick out the ones that are stable with respect 
to changing electron density. Here we try to concentrate 
on a few of these plots and give a thorough discussion of 
them. To this end we calculate the time-dependent, in- 
duced charge density at the absorption-peak frequencies 
in four adjacent unit cells, and plot snapshots thereof in 
Fig. H (b) and (c). The antidots are situated at the inter- 
sections of the thick lines and the ( "-" ) signs mark 
the locations of the charge density maxima (minima). 

The two spectra shown in Fig. ^ (a) calculated for dif- 
ferent electron densities are very similar. At the same 
time, however, the induced charge densities at the differ- 
ent peaks can in general change quite a lot with changing 
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electron density. There are a few exceptions to this, most 
notably the peaks marked (H) and (L) and indicated by 
arrows. As we will see, one can give a clear, semiclassical 
interpretation to these modes. 

Thus, Fig. 1^ (b) shows the charge density correspond- 
ing to the (H) peak in Fig. || (a). One sees a dipole which, 
looking at a sequence of snapshots, rotates around the 
center of each lattice cell (i.e., between four antidots) in 
the direction of cyclotron motion. This mode, which can 
be anticipated on general grounds, emerges in simple tlift^ 
oretical modelsQ and has been detected experimcntallyD'El 

In the low frequency region we find a more complicated 
collective mode [peak (L) in Fig. ^ (a)] depicted in Fig. 
H (c). A dipolar charge distribution is rotating around 
each antidot in the direction of cyclotron motion, and a 
"ring" of three charge density maxima and three minima 
between the antidots is rotating in the opposite direction. 
During each period one sees some small charge transfer 
between the two structures. We interpret this mode as 
a pair of coupled (inter)edge magnetoplasmons with an- 
gular momenta I = +1 and I = —3 around an antidot 
and the center of a cell, respectively. The dynamics of 
this mode is mainly determined by an equilibrium be- 
tween the Lorentz force and restoring electrostatic forces. 
From this follows that the magnetoplasmon propagates 
in opposite directions around a charge density, maximum 
(a cell center) and a minimum (an antidot). Note also 
that this mode is an example of mixing of states with an- 
gular momenta differing by a multiple of four in a square 
lattice. 

Both peaks, (H) and (L), show absorption of light po- 
larized in the direction of the cyclotron resonance in 
agreement mith experiment The existing theoretical 
explanation,!] however, is based on a model with circu- 
larly symmetric unit cells and cannot describe the inter- 
play of modes centered at different places of the unit cell. 
Besides the modes discussed until now the 2DES absorbs 
energy at a number of other frequencies. The correspond- 
ing induced charge distributions are more complex than 
the ones displayed in Figs. || (b) and (c). These exci- 
tations are to a large extent of a quantum-mechanical 
nature, i.e., the result of intersubband and inter-Landau- 
level transitions (see also Ref. In this context, it is 
also clear that our spectra obtained for short-period su- 
perlattices are not completely comparable to the ones 
found experimentally. There are two main reasons for 
this, the potential has a stronger influence on electron 
motion, and we have not treated disorder broadening. 

In conclusion, we have developed a theory that makes 
it possible to study the electronic structure of a 2DES in 
a combined periodic potential and strong magnetic field 
in a detailed fashion treating electron-electron interac- 
tions at the mean-field level. The so calculated energy 
spectra show clear traces of a self-similar structure like 
the Hofstadter butterfly, however, considerably modified 
by screening effects. The dielectric response of the 2DES 
in the FIR regime reveals a rich spectrum of excitations. 
Some of the peaks in these spectra can be interpreted in 



terms of semiclassical collective excitations, while others 
mainly are of a quantum-mechanical origin. 
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FIG. 1. The width of the four lowest Landau bands (in units hcuc) plotted versus inverse magnetic flux. The four principal 
Fourier components of the potential [v^^^iO, ±1) and u'"'*(±l,0)] arc set to 1.43 meV, a = 1000 A, and = 1.2 • 10^^ cm~^. 
The bands are centered around the limiting Landau level values n + 1/2 and get broader as the magnetic field decreases. The 
commensurability phenomena manifest themselves in the intricate splitting of the bands. To underscore the nonmonotonous 
dependence of broadening on the magnetic field, we also display the overall band widths in the inset of the left upper graph. 




FIG. 2. Absorption spectra (panel a) and the oscillating charge-density distributions associated with the peaks indicated 
by arrows (H panel b) and (L ^ panel c), respectively. The upper curve in (a) is offset by 10. The external potential has 
Fourier components v'"'\Q, ±1) = v'"'\±l, 0) = huic = 2.86 meV and v'"'\0, ±2) = i;''''*(±2, 0) = hcu^/A:. For the high-frequency 
mode one sees in (b) a dipole rotating around the center of each lattice cell. For the low-frequency mode one can in (c) observe 
a dipole centered on each of the antidots and a hexagonal pattern between the antidots. One structure of each kind is marked 
with signs. 
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